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Abstract 

Interference phenomena of neutrinos are studied. High energy neutrino 
in T2K near detector and low energy neutrino in KamLAND are possible 
experiments that could show macroscopic interferences of neutrino waves. 
In both experiments interference patterns may give new insights on the 
absolute value of the neutrino mass. 



1 Interference of neutrino waves 

If a wave function tjj u (t,x) is a sum of two different wave functions, the prob- 
ability of observing this particle that is proportional to \ip(t,x)\ 2 has an inter- 
ference term. Interference phenomena of the light, electron, and neutron have 
been studied well and are used for many purposes in wide area. The neutrino 
is a wave and interacts with matters so weakly that it is extremely difficult to 
control the neutrino. It is a challenging task to prepare a detection method for 
the neutrino interference. 

In this article, we show that diffraction-like experiments of neutrinos are 
possible despite its extreme weakness of interactions with matters [TJ [2] and that 
they give new insights on neutrinos. Due to the weakness of the interactions 
it is necessary to analyze whole processes from its production to detection [3]. 
Neutrino becomes a wave packet of a finite coherence length and behaves as a 
wave and particle simultaneously. Its coherence properties are important for 
the observability of the interferences. Neutrino is produced in the scattering 
or decay of a particle by the weak interaction and is in the intermediate state 
between the production and detection processes]!] as is shown in Fig. [TJ 

particle + target — > neutrino — > neutrino + target — > particles. (1) 

For the interference pattern of the neutrino to become observable, a wave at the 
detector should be a superposition of multiple waves of different phases. 
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Fig. 1: The whole processes of the neutrino reaction is shown. The neutrino is 
produced as a decay product of the pion or nucleon and propagates for macro- 
scopic distance. The neutrino is detected finally. 

Distance, L, and energy, E v , of the neutrino propagations of T2K [5 near 
detector and KamLAND [6], which we focus in the present work are 

T2K : L = few x 100 [m].; E v = few x [GeV/c], (2) 
KamLAND : L = 200 [km].; E v = few x [MeV/c]. (3) 

L is the distance between the reactor and detector in KamLAND, and is shorter 
than the distance between pion source and detector in T2K, since pions decay 
in large decay area. These parameters are different each others but two experi- 
ments show similar phenomena of interferences. 

2 Neutrino production and detection :wave packet 
scattering 

The amplitude of the neutrino reactions from its production to detection is 
written as 

(l,a out ,p out \ J dx 1 dx 2 T(H l (x 1 )H i (x2))\a in ,p in ) 7 (4) 

with a weak interaction Hamiltonian Hi(x) and the in-states Q;, n ,/3, n and out- 
states a ou t,/3 ou t and a lepton I. In the standard S-matrix, the in-states and 
out-states are plane waves defined at t = — oo and at t = oo. Hence this 
amplitude is invariant under the translation of the space and time and is useless 
for the space-time dependent informations in scatterings such as probability 
at a finite time interval or finite distance. These translational non-invariant 
quantities are calculated if the states a, /3 are defined using wave packets [3]. 
The wave packet is defined around a certain time and position and has the energy 
and momentum of finite uncertainties and is useful for studying the space time 
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dependent informations. We study the amplitudes of the wave packets and find 
the finite time or length dependences of the amplitudes. 

First the coherence lengths in which the waves keep coherences are evaluated 
from the production and detection processes. Neutrinos are produced by a decay 
of particles which have certain coherence lengths. As was analyzed in Ref.[T], the 
coherence lengths are the sizes of waves that maintain coherence. The mean free 
path, on the other hand, is calculated from the transition probability and the 
density of scatterers as the average distance that a particular wave propagates 
freely. Hence the coherence lengths arc determined by mean free paths. 

(l-l)Decay of pion in fright. 

Pion is produced in matter, in metal for instance, and the mean free path 
was calculated p] as 

Vpion = 10 _1 [m]. (5) 

When these pions are emitted from metal into the vacuum, these waves maintain 
coherence within this size. So the length in which particles maintain coherence 
is determined by the mean free path. These particles are described by the wave 
packets of the finite coherence length in the vacuum. 
(l-2)Decay of Nucleus in solid. 

The other sources of neutrino are unstable nucleus. In reactor, these unstable 
nucleus are bound in atoms, and they are described by wave functions of finite 
sizes. The sizes of nucleus wave functions are estimated from the center of mass 
gravity effect between the nucleus and electrons in the atom Q] as, 

o Nucleus = rx R atom = x 10" 10 « 10~ 13 [m], (6) 

where r is the ratio between the mass of nucleus and electrons and R a tom is 
the size of atom. The size of nucleus wave packet is in the range between the 
nucleus size and atom size. So, the size of pion wave packet in fright is large 
but the size of nucleus wave packet is small. 

Neutrinos interact with particles in the target such as electrons or nucleus 
which have finite coherence lengths. The electron has a large size but the nucleus 
has a small size, and they are defined by wave packets. 

(2-l)Electron wave function in matter. 

Bound electrons in atom are described by atomic wave functions. The wave 
packet of the electron has a size of the electron wave function 

^electron = 10" 10 [m]. (7) 

(2-2) Nuclear wave function in solid 

Nucleus is extremely small and is described by a nucleus wave function. 
From the size of nucleus wave function Eq. ((6]) , the wave packet of nucleus has 
a size 

^Nucleus = 10 _13 [m], (8) 
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which is the same as Eq. ©. Thus the nucleus has a small size. 

Wave packet is a superposition of plane waves with a momentum dependent 
weight. Gaussian and Lorentzian wave packets are studied often. They become 
maxima around the central region and decrease uniformly at large momentum 
region. The former decreases rapidly but the latter decreases slowly. In the cen- 
tral region , Gaussian wave packet is most convenient for practical calculations 
and is applied here also. In the tail region, where the momentum becomes large, 
the particle correlations are determined here from their production processes. 

The tail of the wave packet of pion that is produced in hadron collisions 
depends on the production processes. They are described by pion correlation 
functions 

A p { Xl - x 2 ) = {a\<p{ Xl )\l3){p\<p{x 2 )\a), (9) 
for the initial and final state |a), \/3) and is transformed as 

A„(xi - x 2 ) = / dV^-^Apfo). (10) 



The A(q) of a free relativistic particle, A p ^ ree \q), and of an interacting particle 
with ^(x) interaction, A p nt \q), are 

A(/ree) (?) = 7^4^ V— . ( n ) 

(z7r)^ q z — m z — it 

^ {mt) ^) = TTV i - 1 - T —5((q + qs) 2 ~m 2 ). (12) 

(27r) 4 q z — m z — it q z — m l + xe 

From these correlation functions, A p ^ ree \x) and Ap int ^(a;) are calculated in 



Euclidean metric in the limit A = V x 2 — > as 

A p Sree) (x) = + less singular term, (13) 

A P mt) (*) = ^2+°( A )> : 95^0- (14) 



3 Wave packet evolution 

A wave packet is a superposition of the momentum states and has a finite 
spatial width 7 8 [9 [10 . Although a neutrino wave packet is an intermediate 
state in our formalism, it is instructive and useful to analyze the wave packet 
and its evolution in time here. The wave packet behaves like a particle and has 
the center position of momentum, space, and time. Although it behaves like a 
particle, it also behaves like a wave and has a characteristic phase. So the wave 
packet is not the same as an ordinary particle. 
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Wave packet of the central values of the position, Xq, the momentum, Po, 
and the time, T is decomposed as 

(t,x\P ,X ,T ) = J dp(x\p)(t,p\P ,X ,T ) (15) 

= jV 3 ( — ) * y ^ e -^(p)(*-ro)+ip-(S-^o)-f (P-Po)^ 

where <r is the size of wave packet. The wave packet maintains its shape and 
size during small t — T period, and it has a translational motion. 

(t,x\P ,Xo,T ) =iV 3e -^(^--?o-i7o(t-To)) 2 e ^ ; ( 16 ) 

cf> = E(P )(t - T ) - Po ■ (x-X ), vo = (17) 

The wave is extended around the center x(t) = Xq + Vo(t — To), and is 
classified by the phase 4>. The phase at the center is given by 

4> = E(P )(t - T ) - Po ■ (x-Xo) =(E-^j(t- To) = ^(t - T ). (18) 

This phase is proportional to the mass squared and inversely proportional to 
the energy. Hence its magnitude becomes extremely small. At large t — To, the 
stationary momentum that satisfies 

JL(- lE (p)(t-To)+ip-(x-X )-^(p-Po) 2 ) =0. (19) 

gives a dominant contribution to the above momentum integral. The momentum 
which satisfies this equation is determined by the time and space coordinates, 

(20) 

Pf = m 1 (x-Xq) (21) 

^{t-Tof-{x-XoY 

and the wave function is described as 

(t,x\P ,Xo,T ) (22) 

This wave has the same phase 

0=-g(t-T o ) (24) 
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as before and the magnitude varies with time. The magnitude decreases and the 
wave expands with time in space. The expansion parameters in the longitudinal 
direction, jl, and in the transverse direction, "fx are 



- ' it ' (25) 

The longitudinal expansion is proportional to the mass squared and vanishes in 
the massless case. The transverse expansion is independent from the mass and 
is determined by the initial size of the wave packet. 



4 Probability at finite time interval 

Probability at finite time interval :T2K near detector 

The transition amplitude of the neutrino reaction which includes the pro- 
duction and detection processes are studied by wave packets [IT] . The neutrino 
is treated as a wave in the intermediate states and its size is determined by the 
particle and wave properties in the initial and final states. 

The amplitude T for a T2K neutrino reaction that includes the production 
and propagation processes where the neutrino is in the outgoing state is given 
as 



T = igm^N' / dtd&dp v {/3f\<p(x)\aii) 



N' 



xe! (fi&)t-f^\( p ;) 75 ,( p ;) e <fift)(*-T,)-f.(-i„))-^( P V 



(2tt)5 \Efa). 

3 

(Pf\<p(x)\ai) = ^j 7 e -<sft-)(".)-fe»-(-i)) 

xe -^(-i-^(»-^)) 2 ^, ai . 

In the above equation, \oti) is the initial state and \/3f) is the final state. <p(x) is 
the pion field and (t, x) is the space time coordinates where the weak interaction 
takes place. The pion is expressed by the wave packet of the size ov and is 
produced at the space time coordinate (T„,X V ). The neutrino is detected with 
the wave packet of the size a v and is detected at (T^, Xj,). 

After the tedious calculations, we have integrated probability of observing 
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Fig. 2: The length dependence of the oscillating part in the 1 [GeV] neutrino 
detection probability of the mass 1 [eV/c 2 ] and 0.5 [eV/c 2 ] is shown. 



neutrino at a finite time T as, 

J dpmuondXneutrinodppion ^ ^ |^| 



T e l -Ef(*i-*2) 



N / dhdt 2 

Jo (ti—ti) 



normal 6 



iuJiv(ti — £2) 



[E 



universal 



(0) X g(T, W„) + /normal X g(T, UJ„ + Ul v )] , 



(28) 



Aprob 



(2y/2ir) 3 a v f„ /„ = 5 2 ™ 2 (m 2 - m 2 ), 







^ = — , L = cT, 



(29) 



where the coefficients -^universal is determined from the large momentum con- 
tribution at the tail of the wave packet Eqs. (|14p and /normal is determined from 
the states at the central region of the wave packet. The probability at a finite 
time T has a T-linear term and an oscillating term. 

The oscillating term for the mass 1 [eV/c 2 ] and 0.5 [eV/c 2 ] is given in Fig. 
The probability becomes minimum at around 300 [m] for 1 [eV/c 2 ] and at 
around 600 [m] for 0.5 [eV/c 2 ] and is expected to be observable. 
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KamLAND neutrino:flavour oscillation 

The neutrino production and detection processes are the second order weak 
process at the coordinate x^ and the detection of the coordinate x 2 . The par- 
ticle states oti and j3i in the initial state are defined by the wave packets and 
integration on the variables, x\ and x% are made easily. The amplitude for 
KamLAND neutrino becomes, then, 

t = in' J 'dximfaxi-ri-f 1 j ' d P- 1 i^f e ~ iMXl ' X2) ( 3 °) 

where Xf and Xl^ are the space and time coordinates of the particles at pro- 
duction and detection and T\ and T 2 are the amplitudes where the neutrino is 
produced or detected. The integrand is an exponential function of those that 
have real parts and imaginary parts and the integral is obtained either around 
the minimum of the real part, Gaussian point, k v , or the minimum of the imag- 
inary part, stationary phase point. The wave packet sizes <7j and a® are defined 
from the wave packet sizes of particles. The Gaussian approximation is good at 
small time and the stationary phase approximation is good at large time. 
Gaussian point 

Using integral around Gaussian point, the neutrino in the above amplitude 
is described by a wave function Tq(k u ,X\ — X2) that is described as 

T (t,X 1 -X 2 ) =exv{i<j>{k v ))f{k v ,X 1 -X 2 ), (31) 
f{k v ,X x - X 2 ) = exp [-aJ((J?! - X 2 ff - b L x {{X x - X 2 ) L - v(X° - X°)f ] 

^ 1 ^ 1 

where 4>{k v ) is the phase of the intermediate neutrino. 
(P-detection : large ) 

When <7£ is large, the variable X® — X° is extended in wide area and is 
integrated with the phase e iE ( k ^- E ^)( x i- x 2) . Then, we have E(k v ) = E v and 
the phase of the amplitude becomes 

4> = K-(Xx-X 2 ) 

= E u \(Xi - X 2 )\ + (j>(m), 

9 

771 — t — t 

4>{m) = ^-\( Xl -X 2 \, (32) 

which agrees to the standard formula of the phase of a momentum state. The 
wide wave packet is approximated well with the plane wave and the oscillation 
phase of the wave packet agrees to that of the plane wave[17j. 
( X-detection : Small <tl ) 

In the small ol, the wave function becomes finite only in a narrow region 
around the position, 

Xi = X 2 + v„{X® - X 2 ), (33) 
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and the time x° integration in Eq. (l3"TT) is made by substituting this relation to 
Eq. (f3"Tj) and the phase is rewritten as 

= {-E(k v )— + k v )\X 1 -X 2 \ = -^-\X 1 -X a \, (34) 

Thus the phase of amplitude for the small wave packet is given in Eg. (|3"4"| 
and that for the large wave packet is given in (|32|) . They have different forms. 
In flavour oscillation, the constant phase which does not depend on the neutrino 
mass cancels and the mass dependent phase <j)(m) determines the interference 
term. The mass dependent phase <j>(m) is described by a unified formula, 

2 

TYI — > — > 

^)^-—^-\Xi-X 2 \, (35) 

with a new parameter 7. 7 = 2 is the value of wide wave packet and agrees to 
that of the standard oscillation formula discussed usually and 7 = 1 is the value 
of the narrow wave packet and reveals a phase of the relativistic particle |12) . 
Stationary phase approximation 

In Eq. (|30[) . the momentum integration is made using the stationary phase 
approximation. The stationary momentum, pi°^ , for a = is 

0(0) = m jx-Xy) 

((x°-X2) 2 -(Z-X„y)i 

is proportional to the mass. After the time x° integration is made, the amplitude 
becomes 

T (k,X 1 -X 2 ) = e -*#(* -^)e^(^ (0, -^) 2 \ (37) 



which is the same as the above X-detection Eq. (|34l) . 
Transition from P-detection to X-detection 

Since the oscillation formula for the narrow wave packet is different from the 
wide wave packet, it is worthwhile to know the wave packet size of the boundary. 
In a detection where the wave packet size L w .p is equivalent to de Broglie wave 
length X W) the X-detection should be better than the P-detection. So naively 
the boundary between two detections is expected at 

Lw.p — A lu (38) 

The wave packet sizes for the KamLAND are given in Eq.© and (|H]). Using 
these values we find from a numerical study of the integration of Eq.fJ^H) that 
the transition from the 7 = 2 to 7 = 1 occurs at E v = 2 [McV/c]. This value is 
the end region of KamLAND neutrino detections and the previous analysis of 
KamLAND will not be affected. In other neutrino detectors that use electrons, 
the wave packets are much larger than the de Brogle wave length for the neutrino 
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energy in a few [MeV/c 2 ], hence the detection are regarded as P-detection. Its 
oscillation length is that of the standard formula. 

At the interface region between the X-dctcction and P-detection where the 
amplitude becomes a superposition of those of 7 = 1 and 7 = 2, 



/(mi) = h{mi) + / 2 (m,) 

m. 2 m. 2 

fi(m i )=e % -* t ,Mm i ) = e*** t , 



(39) 



the probability of flavour oscillation of the masses m\ and mi in this region 
behaves as 

2 



D/1W+/2W 



-it 



,'<5m 2 \ „ (bra 2 \ „ / Sm 2 + m 2 \ „ / Sm 2 + m 
= 2 + 2 cos — — t + 2 cos — — t + 2 cos U + 2 cos 
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where 5m 2 = m\ — m\. The oscillation is determined by the mass-squared 
difference Sm 2 and the absolute value of one mass m 2 . For Sm 2 <C m 2 , the 
oscillation period is determined by the absolute value of the mass mi in this 
region, and for Sm 2 w m\ the oscillation becomes complicated functions of the 
time and energy. 



5 Results and summary 

Transition amplitudes of neutrinos where the production and detection pro- 
cesses are fully taken into account are studied. This amplitudes show that the 
macroscopic interference of the neutrino wave is possible. The probability at the 
finite distance show an unusual oscillation that is caused by the interferences 
of neutrino wave produced at different positions. Using T2K near detector, 
this new oscillation might be observable if the mass of the heaviest neutrino is 
around 0.5 — 1 [eV/c 2 ]. 

This amplitude is applied to KamLAND neutrino and generalized oscillation 
formula, 

T (k v , Xl -X 2 ) = e-^-^e-^^-^T, (41) 

is obtained for the neutrino amplitude. The 7 depends on several parameters 
and a transition from 7 = 1 to 7 = 2 is expected at E v = 2 [MeV] . 
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